APS/V2-05/10/04 
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Solar wind turbulence is dominated by Alfvenic fluctuations but the power spectral exponents 
somewhat surprisingly evolve toward the Kolmogorov value of —5/3, that of hydrodynamic turbu- 
lence. We show that at 1AU the turbulence decomposes linearly into two coexistent components 
perpendicular and parallel to the local average magnetic field. The first of these is consistent with 
propagating Alfven wavepackets and shows the scaling expected of Alfvenic turbulence, namely 
Irosnikov- Kraichnan. The second shows Kolmogorov scaling which we also find in the number and 
magnetic energy density, and Poynting flux. 
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The solar wind provides a unique laboratory for 
the study of Magnetohydrodynamic (MHD) turbulence 
with a magnetic Reynolds number estimated to exceed 
10 5 in the solar windpj. In-situ satellite observations 
of bulk plasma parameters strongly suggest the pres- 
ence of turbulence via the statistical properties of their 
fluctuations 0, 0- Quantifying these fluctuations is also 
central to understanding both the transport of solar en- 
ergetic particles and galactic cosmic rays within the he- 
liosphere, and solar wind evolution with implications for 
the mechanisms that accelerate the wind at the corona. 

The observed fluctuations in the solar wind present a 
complex mixture of hydrodynamic and Alfvenic signa- 
tures. Alfvenic fluctuations dominate the power in these 
fluctuations and are observed propagating away from the 
sun implying solar origin (e.g. However the power 

spectra jj, 0, E| suggest an exponent evolving toward 
the Kolmogorov 6] (hereafter K-41) value of ~ 5/3, that 
of hydrodynamic turbulence. This is paradoxical since 
for ideal MHD the turbulent cascade is expected to be 
mediated via Alfven wavepackets suggesting an exponent 
of <~ —2/3, that of Irosnikov and Kraichnan (here- 
after IK). Intervals can be found where different magnetic 
field and velocity components simultaneously exhibit K- 
41 and IK scaling|j| 0, indeed, these phenomenologics 
can be difficult to distinguish in low order moments |lfj|. 
The flow is also observed to be intermittent, this has 
been suggested to account for the 'anomalous' —5/3 scal- 
ing in the power spectra in terms of incompressible MHD, 
rather than hydrodynamic, phenomenology ^lj- Alfvenic 
fluctuations, when isolated by the use of Elsasser vari- 
ables (see e.g.Q), and decomposed by considering dif- 
ferent average magnetic field orientations that occur at 
different times, are found to be multicomponent |l2j| . and 
coupled [l3|. This picture, of an essentially incompress- 
ible, multicomponent Alfvenic turbulence |j, 0] suggests 
that a significant population of Alfvenic fluctuations 
evolve to have wavevectors almost perpendicular to the 
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background magnetic field, leading to a 'fluid- like' (in the 
sense of K-41) phenomenology, and the —5/3 power spec- 
tral slope. However, fluctuations in solar wind den sity 
are not simply proportional to that in magnetic fieldjlj] 
and show nontrivial scaling |2|> Il5| that suggests that the 
turbulence is compressible [lq . The role of compressibil- 
ity is thus an open question. An important corollary is 
that the full behaviour cannot be captured by models 
which describe the observed Alfvenic properties in terms 
of fluctuating coronal fields that have advected passively 
in the expanding solar wind|l7|. 

Here, we will quantify the interplay between K-41 and 
IK phenomenologies in the turbulent solar wind. We 
can discuss the statistical properties of fluctuations in 
some variable of the flow, such as velocity, magnetic 
field, or density, by considering ensemble averages. Fluc- 
tuations in the velocity field can be characterized by 
the difference in some component, or in the magnitude, 
Sv = v(r+L)— v(r) at two points separated by distance L. 
The dependence of Sv upon L is determined in a statis- 
tical sense through the moments < Sv p >, where < ... > 
denotes an ensemble average over r. Statistical theories 
of turbulence then anticipate scaling < 5v p L >^ L^( p \ 

Kolmogorov's 1941 theory for hydrodynamic flows[|J 
essentially follows from dimensional analysis. A fluctu- 
ation Sv arising from a transient structure in the flow 
with characteristic lengthscale L, and timescale T, trans- 
fers kinetic energy Sv 2 implying an energy transfer rate 
€l ~ Sv 2 /T ~ Sv 3 /L. If the statistics of the fluctua- 
tions in the energy transfer rate are independent of L, 
its p moments < e p L >~ 6q where the constant eo is 
the average rate of energy transfer. This gives the K- 
41 scaling < Sv p L >^ L p / 3 . In practice, hydrodynamic 
flows are found to deviate from this simple scaling. This 
intermittency^] is introduced through a lengthscale de- 
pendence of the fluctuations in energy transfer rate so 
that < e p L >~ t p (L/L y(P\ where Lq is some character- 
istic lengthscale and fi(p) is the intermittency correction. 
The scaling for the moments then becomes < 8v p L >^ 
L c(p) with the K-41 exponents C(p) = p/3 - ^{p/3). For 
incompressible MHD turbulence, Alfvenic phenomenol- 
ogy mediates the cascade, introducing an additional char- 
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FIG. 1: Structure function analysis of solar wind density fluc- 
tuations. Inset: structure functions versus differencing inter- 
val r (the traces are offset for clarity). Main plot: scaling 
exponents computed from the raw data (*), and applying 
an upper limit to fluctuation size of 20<t(t)(D), 15ct(t)(o), 
10ct(t)(A) and 5<t(t)(o). 
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FIG. 2: Scaling exponents £(p) versus p for solar wind quan- 
tities, in blue, magnitudes of velocity «(□) and magnetic field 
B(o), in black, number density p(o), number flux pv(A) and 
momentum flux pv 2 (o), magnetic energy density _B 2 (v) and 
Poynting flux vB 2 (M). Note that £(4) « 1 and f(3) w 1 
respectively for these groups of quantities. 



acteristic speed, the Alfven speed va- The above dimen- 
sional argument then gives an energy transfer rate cl ~ 
(8v 2 /T)(Sv/va) ~ Sv^/L, which is just that proposed by 
Iroshnikov and Kraichnanpi , so that < Sv p >^ L^ p ' 
now with £(p) = p/4 — fi(p/4). 

The experimental study of turbulence then centres 
around measurement of the scaling exponents, the C(p)- 
A full description requires the (difficult to determine) 
intermittency correction, the p{p). However, if the sys- 
tem is in a homogeneous steady state, the average en- 
ergy transfer rate is uniform so that < >= eo an d 
//(l) = giving, for K-41 C(3) = 1, and for IK, C(4) = 1, 
independent of the intermittency of the flow. A deter- 
mination of the lower order moments that is sufficiently 
accurate to distinguish these two cases is possible for in- 
situ observations of the solar wind and we present this 
here. These observations are typically time series from a 
single spacecraft so that the ensemble averages that we 
will consider will be over time rather than over space, 
the spatial separation L above being replaced by a time 
interval r- the Taylor hypothesis 0. Consistent with al- 
most all experimental studies of turbulence we consider 
generalized structure functions of a given parameter x: 
S p (t) — (| x(t + t) — x{t) \ p ). Solar wind monitors such 
as the ACE spacecraft spend several-year long periods 
in orbit about the Lagrange point sunward of the earth. 
We analyse 64s averaged plasma parameters from ACE 
for the interval 01/01/1998 - 12/31/2001, this consists of 
~ 1.6 x 10 6 samples and is dominated by slow solar wind. 
Figure 1 shows the procedure for extracting the scaling 
exponents from the data. The inset panel shows the 
structure functions of fluctuations in the density versus 
differencing interval forp =1—4. There is a scaling range 
for timescales of minutes up to a few hours, the timescale 



for large scale coherent structures. This scaling range has 
been shown to extend up to almost three orders of magni- 
tude via Extended Self Similarity (ESS)jl^. The scaling 
exponents, that is, the CO 5 ); where S p (t) ~ t^ p ', are the 
gradients of these scaling regions, and these are shown 
in the main plot. The error bars provide an estimate of 
the uncertainty in the gradients of the fitted lines (lin- 
ear regression error). Finite, experimental data sets in- 
clude a small number of extreme events which have poor 
representation statistically and may obscure the scaling 
properties of the time series. One method H El for ex- 
cluding these rare events is to fix a (large) upper limit 
on the magnitude of fluctuations used in computing the 
structure functions. Importantly, this limit is varied with 
the temporal scale r to account for the growth of range 
with t in the time series. The figure shows the expo- 
nents computed for a range of values for this upper limit 

1/2 

[5,20]cr(r), where a(r) = . We sec that the scaling 
exponents are insensitive to the value of the upper limit 
once a limit is applied and the rare large events are re- 
moved. Above 10cr(r) this process eliminates less than 
1% of the data points. 

We now compare the scaling exponents for different 
scalar quantities in the solar wind flow, for structure 
functions up to p — 4. In Figure 2 we show the scaling 
exponents for fluctuations in the magnitude of velocity v 
and magnetic field B (blue symbols), along with those for 
the number density p, the magnetic energy density B 2 , 
the flux density pv 2 and the Poynting flux in the MHD 
limit vB 2 (black symbols). The scaling exponents for 
these low order structure functions are determined with 
sufficient precision that we can see that in the case of 
the velocity and magnetic field magnitudes, the £(4) ex- 
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FIG. 3: Structure function analysis of <| <5v ■ b | p > . Inset: 
structure functions versus differencing interval (the traces are 
offset for clarity). Main plot: scaling exponents computed 
from the raw data (*), and applying an upper limit to fluctu- 
ation size of 20cr(r)(D), 15a(r)(o), 10ct(t)(A) and 5<t(t)(o). 
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Structure functions S p versus S3 for p 



6 for 



S p =<| <5v ■ b | p >. The traces are offset for clarity. 
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FIG. 4: Scaling exponents ((p) versus p for the structure 
functions of <| 5v ■ b | p > (□) and of the remaining signal 
(o). Note that £(3) m I and ("(4) « I respectively for these 
quantities. 



ponents are clustered about unity, whereas for the other 
quantities, the £(3) exponents are clustered about unity. 
Thus the fluctuations in velocity and magnetic field mag- 
nitude appear to be dominated by Alfvenic processes, in 
the sense of IK phenomenology, whereas the fluctuations 
in the magnetic field energy, plasma and flux densities 
are hydrodynamic- like, in the sense of K-41. Intrigu- 
ingly, these two phenomenologies are coexistent. 

We now introduce an operation that 'filters out' one 
of the relevant physical processes, namely, Alfvenic fluc- 
tuations. We exploit the property that the full non- lin- 
ear MHD equations support large scale Alfvenic fluctua- 



tions which share a basic property of Alfven waves- that 
the velocity perturbation is perpendicular to the back- 
ground magnetic field. In the turbulent flow, the mag- 
netic field also fluctuates, but we can consider a local 
background value by constructing a running average of 
the vector magnetic field over the timescale t'. For each 
interval over which we obtain a difference in velocity 
5~v = v(t + r) — v(t) we also obtain a vector average 
for the magnetic field direction b — B/ | B | from a vec- 
tor sum of all the observed vector values between t and 
t + r', B(*,r') = B(t) + ... + B(t + -r'), with r' centred on 
t. We choose the interval r' = 2t here as the minimum 
(Nyquist) necessary to capture wavelike fluctuations. Ve- 
locity differences <5v which are Alfvenic in character will 
then have the property that the scalar product <5v • b will 
vanish. This condition filters out all those fluctuations 
which generate a velocity displacement perpendicular to 
the local magnetic field, and is thus less restrictive than 
the ElsasserQ variables which select propagating pure 
Alfven waves. In Figure 3 we plot the structure functions 
of the quantity 5v\\ = 5vb, that is, S p =<\ Svh \ p > ver- 
sus r, (inset) and the corresponding scaling exponents, 
the C(p) for the region where S p ~ r^ p ) (main plot) gen- 
erated in the same way as in Figure 1. In Figure 4 we 
compare these exponents with those obtained for the re- 
maining signal, that is, Sv± — y^Sv ■ 5v — (<5v • b) 2 ). 
Remarkably, both these quantities show a clear scaling 
range (which we will verify) with scaling exponents £(3) 
and £(4) close to unity for dv^ and Sv± respectively. 

This result is consistent with the fluctuations in veloc- 
ity being a simple linear superposition that are close to 
(i) parallel to the local background magnetic field and 
sharing hydrodynamic- like scaling, that is K-41, with 
that of the number, magnetic energy and flux densities, 
and (ii) perpendicular to the local background magnetic 
field with the scaling expected of Alfvenic turbulence in 
the sense of IK. These results provide the first unam- 



4 



o 




6 7 

log 10 (S 4 ) 



FIG. 6: Structure functions S p versus 54 for p = 1 — 6 for 

Sp =<l \/(^ v ■ ^ v - ($ v ■ b) 2 ) | p >. The traces are offset for 
clarity. 



biguous ordering of the data with respect to hydrody- 
namic (i.e.Kolmogorov-like), and Alfvenic (i.e. Irosnikov- 
Kraichnan-like) phenomenology in different parameters 
observed at a single location in the solar wind. 

We verify that these quantities indeed show an ex- 
tended scaling region by means of ESS 20] . If the scaling 

is such that the S p ~ Sq^^^ then a plot of S p versus 
S q will reveal the range of the underlying power law de- 
pendence with t. If, as here, one of the are close 
to unity, the ESS plot will in addition provide a better 
estimate of the C(p)- Figures 5 and 6 show S p versus 
53 for Svh and versus S4 for Sv± respectively, and we 
see that there is scaling over several orders of magnitude. 



The slopes of these plots imply exponents that are mul- 
tifractal, that is, quadratic in p, and distinct for the two 
cases. 

The characteristic nature of solar wind turbulence is re- 
vealed to be a coexistence of two signatures. The first of 
these is consistent with Alfvenic turbulence in the sense 
of IK mediated by Alfven wavepackets propagating par- 
allel to the magnetic field. The second, which shows K-41 
scaling, has a compressive component and could couple 
nonlinearly into strongly oblique, almost non propagat- 
ing Alfvenic fluctuations. This clearly elucidates the pre- 
viously proposed multicomponent nature of solar wind 
turbulence and suggests one of two scenarios. The first 
of these is that the turbulent solar wind is comprised of 
two weakly interacting components- one from the pro- 
cess that generates the solar wind at the corona and the 
other that evolves in the high Reynolds number flow. 
Given the evidence for outward propagating Alfven waves 
and evolution toward K-41 scaling, these correspond to 
IK, and K-41 phenomenology respectively and our re- 
sult yields an important insight into the physics of solar 
wind generation. Alternatively, the coexistence of the 
two components is characteristic of the anisotropic na- 
ture of compressible MHD turbulence in the presence of 
a background field, in which case this determination of 
their scaling properties points to an important modifica- 
tion of theories of MHD turbulence. 
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